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We measure the diffusion rate of Cliern-Simons number in tlie (l + l)-dimensional Abelian Higgs model interact- 
ing with a realistic heat bath for temperatures between 1/13 and 2/3 times the sphaleron energy. It is found that 
the measured rate is close to that predicted by the sphaleron approximation at the lower end of the temperature 
range considered but falls at least an order of magnitude short of the sphaleron estimate at the upper end of 
that range. We show numerically that the sphaleron approximation breaks down as soon as the gauge-invariant 
two-point function yields correlation length close to the sphaleron size. 



In the study of thermally induced transitions 
over energy barriers in field theories numerical 
real-time simulations are often the only compu- 
tational tool reliable in a wide range of tem- 
peratures. Determination of the corresponding 
transition rate in a simulation should include its 
proper averaging over the canonical ensemble in 
the phase space of a system in question. If such 
averaging is done by immersing the system in a 
heat bath, the latter should be coupled in a way 
that leaves the bulk dynamical properties of the 
system intact. This is achieved by allowing in- 
teraction with the heat bath only at boundaries. 
The choice of the heat bath is further restricted 
by requiring that its dynamical properties will re- 
semble as closely as possible those of the system 
itself. We are then led to the following realistic 
heat bath (RHB) recipe: linearize the field the- 
ory to be studied in the vicinity of a vacuum and 
use the resulting system as a boundary-coupled 
heat bath [2]. In our earlier work [1] we verified 
that this construction works as intended in the 
simplest case of one-component field in one di- 
mension. Here we report on applying the RHB 
technique to the study of sphaleron transitions in 
the Abelian Higgs model (AHM) in one spatial di- 
mension [3]. Our interest in these processes is due 
to their analogy to baryon-number violating tran- 
sitions in the Standard Model, believed to have 
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played an important role in setting the baryon 
number of the Universe to its present value [4]. In 
particular, AHM possesses an infinite number of 
degenerate vacua separated by finite energy barri- 
ers. The sphaleron solution corresponding to the 
barrier saddle point is known analytically [5-7,9]. 

The lattice AHM Hamiltonian in the temporal 
gauge is in suitably chosen units [5] 
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where tTj and -E'jj+i are canonically conjugate 
momenta of the complex scalar (Higgs) field (fij , 
and the spatial component of the gauge field 
j4jj_|_i, respectively. The lattice spacing is a. The 
Hamiltonian equations of motion following from 
(1) are easy to integrate numerically if a Carte- 
sian representation is used for the Higgs field. For 
this reason we use them for real-time evolution in 
the bulk of an open gauge-Higgs system. We have 
verified that the Gauss' law constraint 

liE:,:+i-E,.^j)=Im{7r,rj) (2) 

is maintained to a high precision for very long 
integration times. 

On the other hand, RHB is most easily con- 
structed in terms of gauge-invariant variables. 
In order to obtain a suitable set of variables 
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the Higgs field is rewritten in polar coordinates: 
(f)j = pjexTp(iaj). We then define = 
"i+i - - aAjj+i, = ^Ejj+i and intro- 

duce canonical momentum Trj" for pj . The corre- 
sponding canonical momentum for Uj is expressed 
in terms of e variables using (2). It is easy to 
see that pj , tTj and bjj+i form two pairs of 

canonical variables whose dynamics is governed 

by 
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obtained from polar-coordinate form of (1) by 
substituting (2). 

As a heat bath we take AHM linearized in 
the vicinity of one of its gauge-equivalent vacua, 
which, as is well known, is a system of two free 
fields: the radial Higgs field g and the gauge 
field e whose masses are \/2 and ^/^, respectively. 
Those are coupled at the boundary site of the 
AHM ((1) or (3)) each to its interacting counter- 
part. I.e. g to p, and e to e. The field equations 
at the boundary are modified compared to those 
in the bulk by introducing terms to account (a) 
for the heat-bath reaction to the boundary mo- 
tion, and (b) for thermal fiuctuations of the heat 
bath. The reader is referred to Ref. [1] for details 
of the heat-bath construction and numerical im- 
plementation. In order to be able to transform 
from polar coordinates, used for RHB, to Carte- 
sian ones, used in the bulk, we need to keep track 
of the angular variable a of the boundary Higgs 
field. This is done with the help of Gauss' law 
which for the linearized system gives 
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The focal point of our study was determin- 
ing the diffusion rate of the Chern-Simons num- 
ber Ncs = (Stt)-!/ A(x)dx. Following Ref. [5], 
we extracted the rate from Acs(Oi ^^e average 
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Figure 1. Temperature dependence of the transi- 
tion rate F. The solid curves correspond to (6), 
with the value of ^ indicated near each curve. 

squared deviation of Nqs for a lag t. For suffi- 
ciently long t a random-walk behavior sets in: 



Acs(0 = TLt, 



(5) 



where F is the transition rate per unit length. 

Our rate measurements are to be compared to 
the prediction of Ref. [6] derived under the as- 
sumption that F can be identified with the rate of 
vacuum-to-vacuum transitions, and that the lat- 
ter involve configurations closely resembling the 
sphaleron: 
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where Egph is the sphaleron energy. As Figure 
1 clearly shows, the values of F at low tempera- 
ture slowly approach those given by the sphaleron 
approximation. The two nearly coincide at the 
lowest temperature considered, /3 = 14. Note 
also that the lattice spacing dependence of the 
rate is absent in the low-temperature regime, in 
agreement with the findings of earlier work [7,5]. 
But the most notable feature of our results is 
their dramatic departure downwards from the 
sphaleron prediction starting at about /3 = 5. 
Our measurements indicate that the rate prac- 
tically does not grow in a range of temperatures 
to the left of /3 = 5 point. This range extends to 
higher and higher temperatures as we reduce the 
lattice spacing. 

It is natural to ask why this strong lagging of 
the measured rate behind the one given by (6) 
begins in the vicinity of /3 = 5. To this end recall 
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Figure 2. Temperature dependence of the corre- 
lation length deduced from (7). The sphaleron 
size is shown by the dashed line for comparison. 

the underlying assumption of (6): Chern-Simons 
number diffusion is dominated by evolution of 
configurations resembling the vacuum into those 
resembling the zero-temperature sphaleron. It is 
clear that the Nqs diffusion can only be described 
in these terms as long as the Higgs field is corre- 
lated on a length scale larger than the sphaleron 
size (2\/2). The corresponding correlation length 
A can be found by measuring a gauge-invariant 
two-point function [8] 



Cji = (?i*(?i;exp 
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Figure 2 shows the values of A obtained by fitting 
(Cj;) to const X exp(— |j — /|/A). As expected, 
the breakdown of the sphaleron approximation 
for the rate occurs as A becomes smaller than 
the sphaleron size. Our measurements of A are 
also consistent with the a dependence of the rate 
at high temperatures. For example, at /3 = 3 
A = 1.47 for ^ = 10, only about 3 times larger 
than the lattice spacing a = 0.5. Hence the field 
strongly fiuctuates at length scales comparable to 
the lattice spacing. 

At this point it is unclear what causes the 
sharp slowdown of the rate growth at the high- 
temperature end of our measurement range. We 
cannot exclude a possibility that at tempera- 
tures in question crossing the barrier between two 
neighboring vacua in the configuration space of 
the model [9] in close vicinity of the sphaleron 



saddle point is still strongly preferred energeti- 
cally, but is already suppressed entropically. It 
is usually assumed [9-12] that at temperatures 
above the sphaleron energy the rate grows like a 
power of the temperature. It could be that what 
we observe at /3 < 5 is a crossover from the expo- 
nential to power-law behavior of the rate. It is not 
clear, however, that at such a crossover the rate 
should stop growing as our measurements seem 
to indicate. We hope to resolve this puzzling sit- 
uation by rate measurements at still higher tem- 
peratures, as well as smaller lattice spacings. 
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